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Abstract 

Isoscalar and isovector particle densities are derived analytically by using the approximation of 
a sharp edged nucleus within the local energy density approach with the proton-neutron asymmetry 
and spin-orbit effects. Equations for the effective nuclear-surface shapes as collective variables are 
derived up to the higher order corrections in the form of the macroscopic boundary conditions. The 
analytical expressions for the isoscalar and isovector tension coefficients of the nuclear surface binding 
energy and the finite-size corrections to the j3 stability line are obtained. 

1 Introduction 

The simple and accurate solution of some problems involving the particle density distributions uses the 
nuclear effective surface (ES) approximation[l, 2, 3, 4]. It exploits the property of saturation of the nuclear 
matter and a narrow diffuse-edge region in finite nuclei. The ES is defined as the location of points of the 
density gradient maximum. The coordinate system related locally to the ES is specified by a distance ^ 
from the given point to the surface and tangent coordinate rj (see Fig. 1). The variational condition of 
the nuclear energy minimum at fixed other integrals of motion within the local energy density theory, in 
particular, the extended Thomas-Fermi (ETF) approach[5, 6] is simplified much in the ^,77 coordinates 
for any deformations by using expansion in small parameter a/ R ^ A^^^^ <^ 1 for heavy enough nuclei 
(a is the diffuse edge thickness of the nucleus, and R its mean curvature radius). The accuracy of the 
ES approximation in the ETF approach was checked [4] by comparing results of the Hartree-Fock (HF) 
and ETF theories based on Skyrme forces[5, 7] without spin-orbit and asymmetry terms. Within the 
ES approximation a rather reasonable agreement of the calculations with the experimental data on a 
mean particle-number dependence on the excitation energies and reduced transition probabilities of the 
low- lying collective states of non-magic nuclei was found[8]. In the present work, we extend the ES 
approach[4] taking into account the spin-orbit and the asymmetry effects in nuclei. 
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Fig. 1: ES and local rj coordinates. The pro- 
file function y = Y(ri) in cylindrical y, z coor- 
dinates is shown schematically by thick solid 
curve [1, 3]; a is the parameter of the nuclear 
diffuse edge. 



2 LOCAL ENERGY DENSITY AND CONSTRAINTS 

We begin with the nuclear energy E within a local energy functional approach[5, 7]: 

+ p+[e+{p+)-e.] + [A + Bp+ + ^{Wp+f + A-{yp-)\ (1) 

where f[p+(r), p_(r)] is the energy density as a function of the isoscalar p+ and isovector p- particle 
densities. It overlaps approximately most of the realistic Skyrmc forces [7]. p± =p„±/3p, X = {N — Z)/A, 
N = J drpn, Z = J drpp, A = N + Z, $ is the Coulomb potential and $ is its average up to a small 
exchange component [1, 2, 10]. As usually, £ of (1) contains the volume, and the surface terms without 
and with the gradient density terms [1, 3, 4], =16 McV is the separation energy per particle and 
^sym=60 MeV is the symmetry energy constant of the nuclear matter. The semiclassical U corrections 
appear through T = h^/18m in the ETF kinetic energy density [5, 6], m is the nucleon mass. In (1), we 
have neglected relatively small isovector (spin-orbit and semiclassical) corrections. The isoscalar surface 
energy density part, independent of the density gradient terms, is determined by the function £_|_(/o_|_) 
which satisfies the saturation condition: 



e+{p) = 0, de+{p)/dp+ = 0, (2) 

where p = 3/47rrQ ~ 0.16 fm""^ is the density of the infinite nuclear matter, ro = R/A^^^ is constant 
independent of A. For the isovector component one has 

^- = ^ (^2 _ ^2 /^^^ - ^ [(1 - P-/P+) $ - (1 - A-) . (3) 

The spin-orbit gradient terms in (1) are defined with a constant: B = —9mW^/16h^, Wo=100 - 130 MeV 
fm^ (sec rcfs. 5,7). 

From the condition of the energy E (1) minimum together with the constraints for the fixed particle 
number A, neutron excess N — Z, and deformation Q of the nucleus[l, 2, 9]: 

A = j dv p+(r), N-Z = Jdr p_(r), Q = J dr p+{r) q{r), (4) 
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one arrives at the variational Lagrange equations: 

- A+ - Aq g = 0, _ A_ = 0. (5) 

Hcrc, A+, A- and Ag arc the corresponding Lagrange multipUers where A+ and A_ are the isoscalar and 
isovector chemical potentials, respectively. 



3 ISOSCALAR AND ISOVECTOR PARTICLE DENSITIES 

In the nuclear volume, up to the second order in — p one gets[l, 3, 4]: 

s+{p+) = -r^{p+-pf , e(w) = ^ = (1 - w)^ w;=^- (6) 



where K is the incompressibility of the infinite nuclear matter. From the Lagrange equations (5) one 



finds for the volume densities p^ ' 



pV"^-P 1 + ^ , p]!'^-plX+-^ . (7) 



^sym 



Small finite-size corrections of the order of a/R ^ aH [H is a mean ES curvature) are determined by the 
surface components of the corresponding chemical potentials: 

a(+) = A+ + 6„-^A'2-|$ + AQ<z(r)~aif~a/i?~^-V3, 

k[-J = X_-b,ymX + ^^^aH. (8) 



For the dimensionless isoscalar density w{x) (6), from the first equation (5), up to the leading order 
m a/R one obtains the ordinary first-order differential equation: 



dM;__ / e{w) _^ _ j ApK 

dx~ '^\lw + ^w^+j' """a' °"Y18&2' 

where (3 = Bp/ A, 7 = T/A'pA. By differentiating equation (9) one finds the boundary condition from the 
definition of the ES: d^w/dx^ = at a; = = 0), 

{wo-pwl-j)e{wo) + wo{wo+pwl+^)(^^j =0, (10) 

V / W—Wq 

together with the condition of the exponentially vanishing the density outside the nucleus: w ~ exp(— a;) = 
exp(— ^/a). Solving the problem (9), (10), one arrives at the solution in the inverse form x{w): 
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With the quadratic approximation (6) for e(u>) one gets the analytical solutions in terms of the algebraic, 
trigonometric and logarithmic functions. For /3 = 7 = it simplifies to w{x) = tanh^ [{x — xo)/2] for 
X < xq = 2arctanh (l/Vs) and zero for x outside the nucleus[4]. In Fig. 2 (left) the influence of the 
scmiclassical correction to w{x) is shown by comparing the 7 = (dashed line) and the "exact" (thin 
solid line) cases. This correction is small everywhere, besides the quantum tail outside the nucleus for 
a;^ 1. Almost the same results one obtains for the SkM* and SLy7 forces [7]. One should also notice a 
rather big effect of the spin-orbit interaction as compared to the simplest analytical solution at /? = 7 = 0. 
We found also a good convergence of the expansion of the e{w) in powers of 1 — w in the density solution 
(11) by comparing the exact numerical function[7] e(w) to its approximate solution (6). The agreement is 
within a precision of the line thickness. Fig. 2 (right) presents a weak sensitivity of the isoscalar solution 
w{x) (11) on the different Skyrme forces with (6) for e{w). 
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Fig. 2: Density w{x) (11) (left) as a function of x = ^/a and its comparison (right) for several Skyrme 
forces[7] (n =5-7,230a and 230b in SLyn) for e{w) (6). 

For the isovector density up to the leading order in a/R, after simple transformations one finds the 
equation and the boundary condition in the form 



dw- 
dw 



= Cs 



€{W) 



1 - 



W_ 



W-{1) = 1, 



W- = 



(12) 



where Csym = d ^/—bsym/^P -A.-. Up to the leading order of the ES approximation in a/R, one obtains 
the analytical solution through the expansion in powers of 1 — iz;. 



W- = w cos [■u(w)] , u{'w) 



1 



1 + 



1 — W 



Csymil + P)+P/'2, 



(13) 



The dependence of the dimensionless isovector density w-{x) (13) on the semiclassical and spin-orbit 
effects versus the corresponding results for the density w{x) (11) are shown in Fig. 3 (left). A weak 
sensitivity of the dimensionless isovector density W- on the choice of the Skyrme forces is seen in Fig. 3 
(right). 
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4 ES EQUATIONS AND LDM BOUNDARY CONDITIONS 



For more exact isoscalar particle density we calculate the main terms of higher order in the parameter 
a/R in the first equation (5). Integrating this equation over the ES in normal-to-smface ^ direction and 
using the equation (9) up to the leading order in a/R, one arrives at the differential equation 



ES 



ep 



dp+ 



3 ( 
TP iP+ 
2p 



P) 



^sym 



2f 

Pi+) = 2a+H. 



iP-f 



(14) 



This equation can be considered with respect to the unknown sought profile shape y = Y{r]) of the 
ES in the cylindrical coordinates y, z with the symmetry axis z (see Fig. 1) through the curvature 
H = {l/Ri + l/i?2)/2 in terms of the main ES curvature radii [1, 3] 



i?i = jCY{r]), i?2 = -£V (d^Y/dv^) , jC=\i + {dY{r])/dr]f] 



1/2 



(15) 




Fig. 3: Isovector density W-{x) (13) compared to the isoscalar one w{x) (11) as a function of .t = ^/a 
within the approximation (6) to e{w) calculated for the same SLy7 forces (left) and W-{x) for several 
Skyrme forces[7] (right). 

Eq. (14) is associated with the macroscopic boundary condition[10, 11, 12] with the isoscalar capilliary 
surface pressure Pi^' which is proportional to the surface tension coefficient cr+: 



0-+ 



dp+ 



(16) 



For more exact isovector particle density similarly one obtains the isovector macroscopic boundary 
condition [12], 



ES 



pi-) 

-'■a t 



: pXA^J = X 



bsym {p-/p-X)-- ($-$) 



(v) 



(17) 
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with the isovector surface pressure 



Pi") = 2a.H, a. « 2A- ' , (18) 



where cr_ is the isovector tension coefficient. 



5 SURFACE ENERGY 



The nuclear energy E = Ey + Eg (1) in the ES approximation is split into the volume, Ey = —by A + 
bsym{N - Z)'^/2A + cZ$/4, and the surface terms: 

E,=aS=(bi+^+bi-^)s/47rrl a = <j++a-, bi^'' = Airrl a±, (19) 

where S is the surface area of the ES. The energy Eg (19) is determined by the sum of the isoscalar 
bi'^^ and isovector bi~^ surface energy constants. These constants are proportional to the same tension 
coefficients a± which appear in (16) and (18) and expressed through the surface pressures in (14) and 
(18), respectively, 

6^ = ^^ J dw^/{w + pw^+^)e{w), (20) 

= 108a-X^^ C Aw ~ {cos[«(w;)] - w siuValw)] u'{w)f , (21) 

Ktq Jq VI + pw 

where q;_ = A-/ A and see (13) for u{w). Simple expressions for constants (20) and (21) in terms of the 
algebraic and trigonometric functions can be easily obtained by calculating explicitly the integrals over 
w in the above equations with e{w) taken from eq.(6). Neglecting relatively small spin-orbit terms and 
semiclassical corrections one finds the approximate relationship between the isovector and the isoscalar 
energy constants, b~ ~ a_ X^bf. 

In Table 1 the analytical, bi'^Jn with the approximated (6) for e{'w), the numerical, bi'^num with the 
exact e(w), and bi'^'^ from ref.7 are shown for all Skyrmc forces. One can see a very good agreement 
between all these calculations, besides of SIII. Modula of the isovector constants for the Lyon Skyrme 
forces SLyn[7] are much larger than for other ones. The precision of the spin-orbit and semiclassical terms 
of the isovcc;tor energy density part (1) is not enough accurate for all considered Skyrme interactions as 
the isovector surface tension cr_ which appears in the surface energy (19) becomes inconsistent with that 
of the capilliary isovector pressure (17), (18). These terms can be improved by fitting to more detailed 
experimental information. 

The (3-stability line is determined by the equivalence of the neutron and proton chemical potentials, 
A_ = An - Ap = 0, and A^^^^ = 2b'^^ H / AnrlpX , according to (8), (17) and (18). With the finite-size 
correction, one obtains X ^ Xq (l — 2bi~\oH/3bsymX§^, where Xq SA^/^e'^ /lOvobsym is the leading 
term[10]. 



6 CONCLUSIONS 

The asymmetry and spin-orbit terms of the energy density within the ETF with Skyrmc forces were 
taken into account analytically by using expansion in a/R ^ 1 of the ES at any deformation. We 
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SIII 


SGII 


SLy230a 


SLy230b 


SLy4 


SLy6 


SLy7 


6+ 

s,an 


17.1 


11.8 


15.3 


17.3 


17.5 


17.5 


17.5 


15.9 


b+ 


18.5 


11.6 


16.5 


18.5 


18.7 


18.7 


18.7 


17.0 


bt 


16.0 


17.0 


14.8 


16.9 


16.7 


18.1 


17.4 


17.0 


^■'^ t/rn 


60.1 


■56.3 


53.7 


63.9 


64.0 


64.0 


63.9 


63.9 




-3.23 


-3.72 


-1.08 


-7.61 


-26.3 


-26.3 


-15.7 


-10.5 



Table 1: The isoscalar (20) and isovector (21) energy surface constants 6^ with Skyrme parameters [7]. 

derived the ordinary first-order equations for the isoscalar and isovector particle densities giving simple 
analytical solutions. When higher order terms are taken into account one gets equations for the moving 
ES in terms of the macroscopic boundary conditions. Expressions for the isoscalar and isovector tension 
coefficients a± in the surface energy were found as those of the macroscopic capilliary pressures of the 
Fermi liquid edge in these boundary conditions. A simple approximate relation of bi~^ to bi'^^ and the 
finite-size ES correction to the /3-stability condition were obtained. Our approach might be helpful as a 
macroscopic part of the nuclear collective dynamical macro-micromodels [11, 14, 15] based on the ETF (or 
LDM) for a "macro" component [8], and the semi- microscopic approach of the Strutinsky shell correction 
method[6, 13] for study of the low-lying collective excitations and fission processes. 
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